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(57) A method of predicting the level of the HbA 1c 
component of glycosylated haemoglobin in a patient's 
blood. A mathematical model is derived which predicts 
the behaviour of the HbA 1c component level relative to 
the blood glucose level using previously measured 
HbA 1c and blood glucose levels. The model is updated 



each time a new HbA 1G level is measured using that 
new measurement and recent new blood glucose level 
measurements. The updated model is then applied to 
predict the HbA 1c level, between measurements of that 
level, using measurements of blood glucose level ob- 
tained since the last HbA 1G measurement. 
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Description 

The present invention relates to the management of diabetes and more particularly to a method and apparatus for 
monitoring the effectiveness of diabetes treatment. 

5 In the treatment of diabetes, a patient is required to regularly check his blood glucose level using a self-testing kit. 

By comparing the result of a self-test with the blood glucose level which he would consider normal, the patient is able 
to estimate the amount of insulin which should be taken in order to bring his blood glucose level back towards that 
normal. Self-testing kits used for this purpose have today become very sophisticated and reliable and provide an 
excellent means for the short term control of diabetes. However diabetic patients can also suffer problems arising from 

10 their condition which only become apparent in the longer term. An individual blood glucose measurement obtained by 
a self-test provides little or no indication of the onset of such long term problems. 

The basic problem which diabetic patients have relates to the transfer of sugar, contained in the blood, across cell 
membranes. This problem in turn makes it difficult for the body to maintain sugar levels in the blood at the correct level. 
Too much blood sugar (e.g. due to the patient injecting too little insulin) and the patient becomes hyperglycaemic whilst 

is too little blood sugar (e.g. due to the patient injecting too much insulin) may cause the patient to become hypoglycaemic. 
In particular, excessive levels of sugar in the blood result in sugar combining with protein to form glycosylated protein. 
Glycosylated protein is substantially insoluble and gives rise to thickening of the walls of veins and arteries, and thick- 
ening of the myelination of nerves. 

One particular form of glycosylated protein is glycosylated haemoglobin. As glycosylated haemoglobin tends to 

20 remain in the blood in the long term, it provides an excellent indication of the level of glycosylated protein in the blood 
and therefore of the effectiveness of the treatment regime which a patient has been following, as well of course as 
indicating how well the patient is following that regime. 

Glycosylated haemoglobin is composed of three components; namely HbA 1A , HbA 1B , and HbA 1c . The HbA 1c level 
in particular is commonly measured by laboratory test in order to provide information on the long term effectiveness 

25 of diabetes treatment. The HbA 1c level reflects the effectiveness of blood glucose treatment over the 6-8 week period 
preceding the HbA 1c measurement. It has been shown that a low level of HbA 1c in a diabetic patient's blood is a good 
indication that the treatment regime is effective and that the risk of secondary problems related to glycosylated hae- 
moglobin is low. The level of HbA-| C in a healthy person's blood is between 4 and 6 % of the total haemoglobin whilst 
in a diabetic person the level may be significantly higher (e.g. greater than 8%). It is generally sought to reduce the 

30 level of HbA 1c in a diabetic patient's blood to between 6 and 7%. 

Due to the often scarce nature of health service resources, and for the sake of convenience and practicality the 
HbA 1c level in a patient's blood is generally tested only every 3 to 4 months. However, given that the HbA 1c level 
provides an indication of the effectiveness of treatment over the previous 6 to 8 weeks, long periods of ineffective 
treatment, and therefore damage to a patient's health, can go undetected with current testing regimes. 

35 The article 'A Theoretical Model to Predict the Behaviour of Glycosylated Hemoglobin Levels' by Kirk W. Beach, 

J. theor. Biol. (1 979) 81 ,547-561 , describes a mathematical model for predicting the level of glycosylated haemoglobin 
from the blood glucose level. This model is however extremely crude and makes use of the simplification that the blood 
glucose level is either constant, changing only by way of a small number of discrete steps, or varying sinusoidally. 
Application of the model to a real patient necessarily involves a great over-simplification of the behaviour of blood 

40 glucose levels. 

It is an object of the present invention to overcome or at least mitigate disadvantages of known diabetes manage- 
ment techniques. 

It is a further object of the present invention to provide a method and apparatus for providing a substantially con- 
tinuous estimate of glycosylated haemoglobin component levels. 
45 According to a first aspect of the present invention there is provided a method of predicting the level of a glyco- 

sylated haemoglobin component in a patient's blood using previously measured blood glucose and glycosylated hae- 
moglobin component levels, the method comprising: 

deriving a mathematical model of the behaviour of the glycosylated haemoglobin component level relative to the 
50 blood glucose level using previously measured levels; 

updating the model when a new glycosylated haemoglobin component level is measured using that new meas- 
urement and recent new blood glucose level measurements; and 

applying the mathematical model to predict the glycosylated haemoglobin component level, between measure- 
ments of that level, using measurements of blood glucose level obtained since the last glycosylated haemoglobin 
55 component measurement. 

Typically ; blood glucose level measurements are made at a considerably higher frequency than glycosylated hae- 
moglobin component measurements. The method of the present invention may therefore be used to predict the current 
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glycosylated haemoglobin component level in a patient's blood using blood glucose level measurements obtained 
since the last glycosylated haemoglobin component level measurement. As the model is updated each time a new 
HbA 1c measurement is made, the model is capable of tracking changes in the physiology of the patient which cause 
the behaviour of the HbA 1c level to change with respect to the blood glucose level. Changes in the blood glucose 
5 measurement pattern, i.e. the times at which the patient makes blood glucose measurements, can also be accounted 
for. 

Preferably, the mathematical model is a parametric model or a semi-parametric model, where the model is defined 
by one or more model coefficients and a model equation which relate blood glucose level to the glycosylated haemo- 
globin component level. More preferably, the model equation relates the glycosylated haemoglobin component level 
10 to one or more parameters which describe, at least in part, the behaviour (e.g. distribution) of the blood glucose level 
over a preceding ; relatively short, time interval. 

The model equation may be a linear equation in which case said model coefficients are the linear coefficients of 
the equation. The linear equation is of the form: 

y= P 1 Aj 1 + p 2 h 2 +...+ p q h q + c 

where y is the predicted glycosylated haemoglobin level, p are the linear model coefficients, h are the parameters 
which describe blood glucose level behaviour, and cis a constant. 
20 Preferably, the behaviour of the blood glucose level over said short time intervals may be described using one or 

more gaussian functions which model the distribution of blood glucose level measurements. Said one or more param- 
eters (h) may be chosen from the mean, variance, and amplitude of the gaussian function(s) or may be derived there- 
from. 

In the case of a parametric or semi-parametric model, the model may be updated following each glycosylated 
25 haemoglobin component level measurement by recalculating said model coefficients (p). In an alternative embodiment 
of the present invention, the coefficients of the parametric model are adapted following each new glycosylated hae- 
moglobin level measurement using an adaptive algorithm. One suitable adaptive algorithm is Widrows algorithm. Such 
adaptive algorithms are arranged to reduce the error between the predicted glycosylated haemoglobin level and the 
measured glycosylated haemoglobin level. 
30 The glycosylated haemoglobin component predicted using the method of the above first aspect of the present 

invention is one of HbA 1A , HbA 1B , and HbA 1c . Preferably however, the predicted component is HbA 1c . 
The method of the above first aspect of the present invention may comprise: 

transmitting measured glucose levels via a wireless data transmission link from a remote station, available to the 
35 patient, to a central data processing station; 

carrying out said steps of deriving the mathematical model, updating the model, and predicting the glycosylated 
haemoglobin component at the central processing station using the transmitted measurements and previous meas- 
urement data stored at the central station; and 

transmitting predicted glycosylated haemoglobin component levels back to the remote station via the wireless data 
40 transmission link. 

Preferably, the remote station is a mobile telephone or a two-way pager whereby the wireless data transmission 
link is provided by a mobile telephone network. For example, the telephone network may be a GMS network and the 
data may be transmitted by the short message service (SMS). It is to be understood that the term 'mobile telephone' 
45 as used here refers to any portable device which utilises wireless telephonic communication including conventional 
cellular telephones and combined cellular telephone/personal data assistant (PDA) devices. 

Alternatively, the method may comprise carrying out all of the steps of the method of the above first aspect of the 
invention in a portable monitoring device. 

It will be appreciated that the method of the present invention may be applied to blood taken from a human or 
50 animal patient. 

According to a second aspect of the present invention there is provided a method of predicting the level of a 
glycosylated haemoglobin component in a patient's blood, the method comprising the steps of: 

transmitting a blood glucose measurement from a remote station to a central processing station; 
55 predicting at the central station a glycosylated haemoglobin component level for the patient's blood using said 

transmitted measurement and a mathematical model of the glycosylated haemoglobin component level relative to 
blood glucose level; and 

transmitting the predicted glycosylated haemoglobin component level from the central station to the remote station. 
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According to a third aspect of the present invention there is provided a system for predicting the level of a glyco- 
sylated haemoglobin component in a patient's blood using previously measured blood glucose and glycosylated hae- 
moglobin component levels, the system comprising: 

5 means for deriving a mathematical model of the behaviour of the glycosylated haemoglobin component level rel- 

ative to the blood glucose level using previously measured levels, and for updating the model when a new glyco- 
sylated haemoglobin component level is measured using that new measurement and recent new blood glucose 
level measurements; and 

memory means for storing said model and/or the updated model and measured blood glucose and glycosylated 
10 haemoglobin component levels; 

means for applying the mathematical model to predict the glycosylated haemoglobin component level, between 
measurements of that level, using measurements of blood glucose level obtained since the last glycosylated hae- 
moglobin component measurement. 

75 in a preferred embodiment of the system of the above third aspect of the present invention, said means for deriving 

and applying and said memory means are provided by a central computer, the system further comprising a mobile 
telephone or two-way pager for conveying measurement data to the central computer via a wireless data transmission 
link. The system may be further arranged to convey predicted glycosylated haemoglobin component levels from the 
central station to the mobile telephone or two-way pager. 

20 For a better understanding of the present invention and in order to show how the same may be carried into effect 

reference will now be made, by way of example, to the accompanying drawings, in which: 

Figure 1 shows HbA 1G level measurements together with blood glucose level measurements obtained from a 
diabetic patient over a period of 640 days: 

25 

Figures 2Ato 2D show respective sets of blood glucose level measurements, extracted from the measurements 
shown in Figure 1, for periods preceding four HbA 1c level measurements; 

Figures 3A to 3D show respectively predicted blood glucose level distributions for each of the measurement win- 
30 dows of Figures 2A to 2D; 

Figure 4 is a flow chart illustrating a method used to predict HbA 1c levels; 

Figure 5 shows measured HbA 1c levels of Figure 1 together with predicted HbA 1G levels; and 

35 

Figure 6 illustrates schematically apparatus for implementing the method of Figure 5. 

The method described below for predicting HbA 1c levels in the blood of a patient is illustrated using actual blood 
glucose and HbA 1c measurements obtained from a human patient. Figure 1 shows the measured blood glucose and 

40 HbA 1c levels plotted against time (measured in days) over a period of 640 days. As discussed above, HbA 1c meas- 
urements are typically made at relatively infrequent intervals (every 12 to 16 weeks) whilst blood glucose level meas- 
urements are made much more frequently (typically two to five times a day). 

Let Y be the set of N available HbA 1c measurements y n shown in Figure 1, where Y= {y^..-y^} and each meas- 
urement corresponds to a measurement time t = T n , (n = 1 to N). Similarly, let Xbe the set of M available blood glucose 

45 measurements x m also shown in Figure 1 , where X- [x v . .x M } and each measurement corresponds to a measurement 
time t = t m , (m = 1 to M). The set of blood glucose measurements Xis divided into sub-sets X n , each sub-set being 
associated with the corresponding HbA 1c measurement V^and containing those values obtained in the period t> T n 
-65 days to 7^. It is noted that when the interval between consecutive HbA 1c is less than 65 days, adjacent sub-sets 
X n will overlap. For the data of Figure 1 , with twelve HbA 1c measurements, twelve sub-sets of blood glucose meas- 

50 urements are obtained, the first four of these being shown in Figures 2A to 2D for the purpose of illustration. 

A suitable semi-parametric model is applied to each sub-set of blood glucose level measurements X n to model 
the distribution of the measurements within the set (i.e. to model the probability density function). The model is chosen 
to have a relatively small number of parameters. In the present example, the distribution of each sub-set X n is modeled 
using the sum of a pair of gaussian functions: 

55 

P(l) = kG(l\\L l ,a ] ) + (1 - k)G(l\n 2 ,a \) (1 ) 
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where k is a coefficient and 0 < k < 1 , / is the blood glucose level, and the gaussian function G has the form: 



The means (ill-, ,0^) and variances (a^,a|) of the gaussian functions, as well as the coefficient k, are obtained using the 
expectation-maximization algorithm (see Dempster A.P., Laird N.M., and Rubin, D.B. (1977); 'Maximum Likelihood 

10 from Incomplete Data via the EM Algorithm'; J. Royal Statistical Soc; B 39 (1); 1-38). The modeled distribution for 
each of the sub-sets of blood glucose measurements are shown in Figures 3A to 3D. 

Each model can be represented by a parameter vector © n {u^ \i n2 , o£ ., o^, k}. A suitable combination of pa- 
rameters is selected from the parameter vector © n and is used to form a model vector h n at time T n . In the present 
example the parameters m, \a 2 > and a constant 1 are selected, corresponding respectively to , h n2 , h n3 and 

is of the model vector. The model vectors h n and the measured HbA 1G levels y n are used to form the rows of a model 
specification matrix Hand an observation vector z respectively and the linear problem z= Hp constructed, where p is 
a regression coefficient vector. This problem can be more fully written as the matrix equation: 
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(3) 



An estimate of the regression coefficient vector p can be determined from: 

p=H 1 z (4) 

where H" 1 is the inverse of the matrix H (or pseudo inverse in the case that H is not square or is singular) and can be 
obtained using a singular value decomposition technique (Press, W.H. Teukolsky, S.A. Vetterling, W.T Flannery, B.P 
1992, 'Numerical Recipes in C: The Art of Scientific Computing', 2nd ed. Cambridge University Press). The estimated 
regression coefficient vector p can then be used to predict a future HbA 1c level y N+1 from the model vector derived 
from the sub-set of blood glucose measurements X N+1 , i.e.: 

3Vi = P^ Vn .1 + Pz h N^ ,2 + Ps h N^ ,s + Pa Vh ,4- ( 5 ) 

Figure 4 is a flow diagram outlining the main steps of this method of predicting HbA 1c levels. 

As has been explained above, changes in the physiology of a patient are likely to change the behaviour of the 
HbA 1c level with respect to blood glucose level. It is therefore important that the coefficients of the regression coefficient 
vector p be updated regularly. This can be achieved by recalculating the vector p every time a new HbA 1c measurement 
is made. Typically, for a patient with no previous records, the regression coefficient vector is first derived using data 
obtainedf rom a number of other patients. Each time a new HbA-| c measurement is made, the new HbA 1c measurement, 
together with the blood glucose measurements obtained since the last HbA 1c measurement, is used in the recalculation 
of the regression coefficient vector. Assuming that the value of N is maintained constant by removing the earliest 
obtained measurements in turn, after N HbA 1c measurements have been obtained from the subject patient, the esti- 
mated regression coefficient vector will be derived solely from measurements made on the subject patient. 

In order to improve the accuracy with which the HbA 1c level can be predicted, and more particularly to avoid biasing 
of the result by rogue measurements, an estimated regression coefficient vector p is in fact obtained for N different 
data sets, omitting in turn each of the elements y, to y N o\ the above matrix equation together with the corresponding 
row of the matrix H. The mean of the estimated regression coefficient vectors is then obtained after removing the top 
and bottom 10% of each of the vector components. For the data set shown in Figure 1 , the final estimated regression 
coefficient vector is the trimmed mean of twelve estimated vectors p. The coefficients of this final regression coefficient 
vector can then be used to predict the current HbA 1c level using equation (5) above. Figure 5 shows HbA 1c levels 
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predicted using a model obtained with this method together with the measured HbA 1c levels used to formulate the 
model (and as shown in Figure 1). It can be seen that the predicted values closely follow the measured values. 

Further improvements in prediction accuracy may be achieved by increasing the number of parameters p, by 
increasing the size of the model vector h n , e.g. to additionally include 8 2 . Alternatively, a set of HbA 1c estimates may 

s be made using a corresponding set of parametric models, each model being obtained using a different set of param- 
eters, e.g. ill-,, u^o^/cor u^, u^c^/c. Accuracy may also be increased by appropriately transforming (e.g. log, expo- 
nential) one or more of the parameters. 

It is well known that blood glucose levels in a patient fluctuate periodically according to the time of day and more 
particularly according to whether or not the patient has just eaten. It may therefore be more accurate to separately 

10 model blood glucose levels for different times of day. This technique requires the separation of all blood glucose levels 
into separate sets. There might, for example, be nine sets, i.e. measurements made before and after breakfast, lunch, 
dinner, and evening snack. For each set j, a set of parameter vectors 0^ is obtained, where n = 1 — /V. Again, for each 
value of n, a set of suitable parameters are selected to form the model vector h n . The estimated regression coefficient 
vector p can then be obtained as described above. 

75 Given that the calculation of the regression coefficient vector may involve a considerable amount of stored data 

(e.g. up to two years of measurement data), and hence present a considerable degree of computational complexity it 
is desirable to store the measurement data and perform computations at a central data processing unit available to 
many patients. In addition, this arrangement improves data security and, importantly, allows the algorithm by which 
HbA 1G levels are predicted to be easily and quickly updated, should that prove necessary. 

20 Figure 6 illustrates a diabetic management system where each patient in the system is provided with a portable 

electronic blood glucose measurement unit 1 of known type. This unit is modified for coupling measured blood glucose 
levels via an interface 2, for example a notebook computer or PC, to a mobile telephone 3. The measurement data is 
then transmitted via the mobile telephone network to a central data processing station 4 of the type described above. 
Based on recently received blood glucose measurements, and a precalculated regression coefficient vector, the central 

25 processing station 4 estimates the patient's HbA 1c level. The regression coefficient vector can be updated each time 
a HbA-| c measurement is made on the patient either using data entered directly into the central data processing station 
or transmitted to it from, for example, a doctor's surgery. The predicted HbA 1c level can be almost immediately trans- 
mitted back to the patient's mobile phone 3, via the telephone network, where it can be displayed on the phone 3 or 
the interface 2. By integrating the functions of the mobile phone 3 ; the interface 2, and the measurement unit 1, the 

30 patient can be provided with a truly portable HbA 1c level monitoring system. 

It will be appreciated by the skilled person that other modifications may be made to the above described embod- 
iment without departing from the scope of the present invention. For example, the probability density function of each 
set of blood glucose levels X n may be modeled using one or more gamma functions as an alternative to gaussian 
functions. 

35 

Claims 

1. A method of predicting the level of a glycosylated haemoglobin component in a patient's blood using previously 
40 measured blood glucose and glycosylated haemoglobin component levels, the method comprising: 

deriving a mathematical model of the behaviour of the glycosylated haemoglobin component level relative to 
the blood glucose level using previously measured levels; 

updating the model when a new glycosylated haemoglobin component level is measured using that new meas- 
45 urement and recent new blood glucose level measurements; and 

applying the mathematical model to predict the glycosylated haemoglobin component level, between meas- 
urements of that level, using measurements of blood glucose level obtained since the last glycosylated hae- 
moglobin component measurement. 

50 2. A method according to claim 1 , wherein the mathematical model is a parametric model or a semi-parametric model 
and is defined by one or more model coefficients and a model equation which relate blood glucose level to the 
glycosylated haemoglobin component level. 

3. A method according to claim 2, wherein the model equation relates the glycosylated haemoglobin component level 
55 to one or more parameters which describe behaviour of the blood glucose level over a preceding, relatively short, 

time interval. 

4. A method according to claim 3, wherein the model equation is a linear equation having the form: 
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y=p 1 /7 1 +p 2 h 2 +...+ p q h q + c 

where y is the predicted glycosylated haemoglobin level, pare the linear model coefficients, h are the parameters 
5 which describe blood glucose level behaviour, and c is a constant. 

5. A method according to claim 3 or 4, wherein the behaviour of the blood glucose level over said short time intervals 
is described using one or more parametric functions which model the distribution of blood glucose level measure- 
ments and said one or more parameters (h) are selected from or derived from the parameters defining these 

10 functions. 

6. A method according to claim 5, wherein the or each said function is a gaussian function and said one or more 
parameters (h) are chosen from or are derived from the mean, variance, and amplitude of the gaussian function(s). 

15 7. A method according to any one of claims 2 to 6, wherein the model is updated following a glycosylated haemoglobin 
component level measurement by recalculating said model coefficients (p). 

8. A method according to any one of claims 2 to 6, wherein the coefficients of the parametric model are adapted 
following a new glycosylated haemoglobin level measurement using an adaptive algorithm. 

20 

9. A method according to any one of the preceding claims, wherein the glycosylated haemoglobin component pre- 
dicted is HbA 1c . 

10. A method according to any one of the preceding claims and comprising: 

25 

transmitting measured glucose levels via a wireless data transmission link from a remote station, available to 
the patient, to a central data processing station; 

carrying out said steps of deriving the mathematical model, updating the model, and predicting the glycosylated 
haemoglobin component at the central processing station using the transmitted measurements and previous 
30 measurement data stored at the central station; and 

transmitting predicted glycosylated haemoglobin component levels back to the remote station via the wireless 
data transmission link. 

11. A method according to claim 10, wherein the remote station is a mobile telephone or a two-way pager whereby 
35 the wireless data transmission link is provided by a mobile telephone network. 

12. A method of predicting the level of a glycosylated haemoglobin component in a patient's blood, the method com- 
prising the steps of: 

40 transmitting a blood glucose measurement from a remote station to a central processing station; 

predicting at the central station a glycosylated haemoglobin component level for the patient's blood using said 
transmitted measurement and a mathematical model of the glycosylated haemoglobin component level relative 
to blood glucose level; and 

transmitting the predicted glycosylated haemoglobin component level from the central station to the remote 
45 station. 

13. A system for predicting the level of a glycosylated haemoglobin component in a patient's blood using previously 
measured blood glucose and glycosylated haemoglobin component levels, the system comprising: 

50 means (4) for deriving a mathematical model of the behaviour of the glycosylated haemoglobin component 

level relative to the blood glucose level using previously measured levels, and for updating the model when 
a new glycosylated haemoglobin component level is measured using that new measurement and recent new 
blood glucose level measurements; 

memory means (4) for storing said model and/or the updated model and measured blood glucose and glyco- 
ls sylated haemoglobin component levels; and 

means (4) for applying the mathematical model to predict the glycosylated haemoglobin component level, 
between measurements of that level, using measurements of blood glucose level obtained since the last gly- 
cosylated haemoglobin component measurement. 
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1 4. A system according to claim 1 3, wherein said means for deriving and applying and said memory means are provided 
by a central computer (4), the system further comprising a mobile telephone or two-way pager (3) for conveying 
measurement data to the central computer (4) via a wireless data transmission link. 

15. A system according to claim 14 and comprising a blood glucose measurement unit (1) coupled to the mobile 
telephone or two-way pager for transmitting blood glucose level measurements to the mobile telephone or two- 
way pager. 

16. A system according to claim 14 or 15 and which is arranged to convey predicted glycosylated haemoglobin com- 
ponent levels from the central station (4) to the mobile telephone or two-way pager (3). 
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